Persistent Rabi oscillations probed via low-frequency noise correlation 
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The qubit Rabi oscillations are known to be non-decaying (though with a fluctuating phase) if 
the qubit is continuously monitored in the weak-coupling regime. In this paper we propose an 
experiment to demonstrate these persistent Rabi oscillations via low-frequency noise correlation. 
The idea is to measure a qubit by two detectors, biased stroboscopically at the Rabi frequency. The 
low-frequency noise depends on the relative phase between the two combs of biasing pulses, with 
a strong increase of telegraph noise in both detectors for the in-phase or anti-phase combs. This 
happens because of self-synchronization between the persistent Rabi oscillations and measurement 
pulses. Almost perfect correlation of the noise in the two detectors for the in-phase regime and almost 
perfect anticorrelation for the anti-phase regime indicates a presence of synchronized persistent Rabi 
oscillations. The experiment can be realized with semiconductor or superconductor qubits. 



The puzzle of the quantum state collapse due to mea- 
surement [l| is becoming accessible for the experimen- 
tal study in solid state systems. Three experiments on 
non-projective collapse 2, 3] have been recently realized 
with superconducting qubits. These experiments (as well 
as the experiment proposed in the present paper) touch 
upon the most intriguing property of quantum measure- 
ment: the presence of a "spooky" quantum back-action 
[3| , which changes the system to agree with the observa- 
tion, and cannot be explained in a realistic way, i.e. by 
using the Schrodinger equation. 

The quantum coherent (Rabi) oscillations in solid-state 
qubits are usually measured in an ensemble-averaged way 
H 0] and decay within a short timescale, even though it 
can be much longer than the oscillation period. However, 
for a continuous weak measurement of a single qubit, the 
Rabi oscillations are non-decaying and can in principle be 
monitored in real time, as follows, e.g., from the quan- 
tum Bayesian formalism [7] , which is generally similar to 
the formalism of quantum trajectories Q. Persistence 
of the Rabi oscillations in this case is due to the quan- 
tum back-action, which tends to increase the amplitude 
of the oscillations to 100%, thus competing against de- 
coherence. The persistent Rabi oscillations lead to the 
spectral peak of the detector signal at the Rabi frequency 
9, [l(| , which has been recently observed experimentally 
3| (see also [IH). In the present paper we will discuss 
another way of demonstrating these oscillations. 

For definiteness let us discuss a "charge" qubit made 
of a double quantum dot (DQD) populated by a single 
electron, the location of which is continuously measured 
by a nearby quantum point contact (QPC). Analogous 
setups can be realized with spin-based or superconduct- 
ing qubits. The continuous qubit evolution due to the 
quantum "informational" back-action can in principle be 
verified in a direct experiment Q; however, it would re- 
quire high-bandwidth recording of the detector signal (in- 
cluding shot noise) and fast qubit manipulation, that is 
still a big challenge for a real experiment. A simpler 
way to study the back-action is to measure the qubit 
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FIG. 1: Analyzed system: a double-quantum-dot qubit mea- 
sured by two QPC detectors, which are biased by combs of 
short voltage pulses with frequency fl coinciding with the 
Rabi frequency Qr. 



by two detectors [13], so that the first (short) measure- 
ment causes a partial collapse of the qubit state, and then 
after a controllable qubit evolution the second detector 
measures the resulting state. Performing the experiment 
many times and selecting a certain result of the first mea- 
surement, it is possible to find the back-action evolution 
experimentally and compare it with the theory. The same 
idea with a different post-processing (selecting the result 
of the second measurement) has been recently used to 
propose an experiment on weak values [ljj . 

The proposal of Ref. 12| still suffers from very weak 
signals produced by two single-shot measurements. An 



obvious way to increase the signal is to average it over 
a long comb of measurement pulses, but in this case the 
selection of a certain result becomes impossible. For- 
tunately, there is a way to overcome this dilemm a by 
combining the ideas of two-detector measurement [12j . 
persistent Rabi oscillations @, @, and stroboscopic 
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quantum non-demolition (QND) measurement 

We propose to use the following setup (Fig. 1). A qubit 
is measured by two detectors, which are biased with two 
combs of short pulses, so that between the pulses the 
qubit undergoes free evolution due to the Rabi oscilla- 
tions. The frequency of pulses Q coincides with the Rabi 
frequency Qr (one pulse per period in each detector) to 
realize the QND regime [15]. When the two combs are 
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not shifted in time relative to each other, the phase of 
the Rabi oscillations is attracted to one of the two stable 
values, corresponding to the qubit being in either local- 
ized state |1) or \2) at the time of measurement. This 
happens because of the usual collapse in the QND frame 
l4 Il5f and is somewhat similar to what happens with a 
parametrically excited swing. However, because of vari- 
ous imperfections (extra decoherence, etc.) there will be 
switching between the two stable regimes, which leads to 
the telegraph noise in the currents through both detec- 
tors. Even if the experimental measurement bandwidth is 
not wide enough to resolve the switching events (which is 
likely for a present-day experiment), the telegraph noise 
is measurable at low frequency via its spectral density, 
which greatly exceeds the shot noise. 

The telegraph noise originates from the presence of two 
quasi-stable regimes of oscillations because of the QND 
measurement. However, for a significant phase difference 
(p between the two combs of the measurement pulses, the 
measurement is no longer QND, and the telegraph noise 
disappears, so that the low frequency noise reduces to a 
much smaller shot noise. For the QND regime is 

restored again and the telegraph noise reappears. So the 
telegraph noise is maximum at both ip — and <p = tt; 
however, for (p = it is almost fully correlated [16| in the 
two detectors (because in the stable regime both detec- 
tors see the same qubit state: either |1) or |2)), while for 
p = n the detector noises are almost fully anticorrclated 
(because when one detector measures state |1), the other 
detector measures |2) half a period later). Experimen- 
tal observation of such noise correlation/anticorrelation 
would show the presence of persistent Rabi oscillations. 

For the quantitative analysis let us use the quantum 
Bayesian formalism 0] and represent the QPC currents 
I a {t) and Ib(t) as (n — a or n = b) 



In(t) = [I ,n + (A/ n /2) Z(t)] f n (t) + £ n (t) y/\fn(t)\, (1) 

where z = Tr(a z p) is the measured z-component of the 
qubit Bloch vector, f n (t) is the dimensionless shape of 
the comb of measurement pulses for nth detector (/ is 
proportional to the QPC voltage), I , n and AI n are the 
detector average current and response for /„ = 1, and 
£ n (i) is the shot noise with the (one-sided) spectral den- 
sity S n . We assume zero temperature. For rectangular 
measurement pulses we use f a (t) = 1 if \t — IT\ < 8t a /2 
and f a (t) — otherwise, where St a is the pulse duration, 
T = 27r/57 is the comb period, and I is an integer. Simi- 
larly, f b (t) = 1 if \t-lT-tp/n\ < St b /2; this describes the 
comb of pulses of duration dtb with the same frequency 
fi but shifted by the phase ip. The qubit Hamiltonian 
-ffqb = (f2fl/2)cr x describes Rabi oscillations with the 
frequency Qr about the x-axis, and we also assume pure 
dephasing of the qubit (not related to the measurement) 
with rate 7. In order to use the Markovian approxima- 
tion, we assume sufficiently high QPC voltages during 
the pulses Q, and also assume |A7„| <C | -fo.nl- 



The spectral densities of two detector noises S aa (ui) 
and Sbb(u) as well as the cross-correlation noise S a b(w) 
can be found via the Fourier transform S nm (uj) — 
2 JToo K nm (T,t ) e~ lulT dT, where K nm (T,t ) = (I m (t + 
T )In{to)) — (I m {to + t)) (I n (to)) is the correlation func- 
tion and the averaging K nm (T,to) is over to within one 
period T; the averaging is necessary because of periodic 
time dependence f n (t) in Eq. (p}. To find K nm {T,to) 
for t > we use the method developed in [loj ]. which 
essentially follows from the quantum regression theo- 
rem [l7j |. Expressing this correlator as K nm (r,to) = 
(AI n AI m /4:)f m (t + T)f n {to)K zz (T,t Q ), we calculate the 
operator-symmetrized zz-correlator K zz (r,t ) as 

K"(T,h) = E (^(*o)|i><ik z |i)Tr[a z pl i >' to (to + r)], 

i=l,2 

(2) 

where p(t ) is the qubit density matrix at time t 
[(i\p(t )\i) — 1/2 because of the symmetry], while 
p\ l )' to (t + r) is the density matrix at time t + r for 
the qubit starting at time to in the state which 
is an eigenstate of a z . Here we should assume the 
ensemble-averaged qubit evolution, for which the mea- 
surement process is represented by pure dephasing with 
rates /„(t)7™ cas , where 7 ™ cas = (A/„) 2 /4S', i . Notice 
that the same technique works also for two detectors 
measuring different observables: one of them defines the 
starting eigenstates, and the other one enters into the 
trace. While so far t > was assumed, we find the 
correlation function at negative t using the symmetry 
K nm {T,to) = K mn {—T,to + t), and also add the shot 
noise contribution 5 nm S(r)f n (to)S n /2 near r = 0. Be- 
sides this quantum method, we also use below the lan- 
guage of a simple semiclassical analysis. 

Let us first assume short measurement pulses, St a ,b <C 
T = 2ir/Q, exactly matched frequency, f2 = f2#, almost 
no phase shift, \<p\ <C 1, and almost negligible extra de- 
phasing, 7 <gc flu. Then we have usual stroboscopic 
QND measurement 14|, [l5( insensitive to the free evolu- 
tion, and therefore the qubit state eventually collapses 
to either |1) or \2) at the measurement moments. This 
obviously leads to non-decaying Rabi oscillations with 
100% amplitude, which are phase-locked with the mea- 
surement combs (though with a random choice of the 
stable phase). The synchronization happens within the 
QND collapse timescale t ~ t co \ = T/(M a + Mb), where 
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8t n (AI n ) 2 / AS n . We assume M n < 1, 



so that t co i >• T, while 7™£ as T are not necessarily small. 
In the ideal QND case the phase of the Rabi oscilla- 
tions is fixed forever after this gradual collapse; however, 
in a realistic case there will be switching between the 
two regimes (state |1) or |2) at the measurement mo- 
ments) with a calculated below rate T$, the same for 
both switching processes because of the symmetry. If 
we assume rare switching, rsi co i <C 1, then the detec- 
tor current /„ averaged over a coarse graining timescale 
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longer than T and i co i, switches between the two lev- 
els, (7o,n i AI n /2)(5t n /T), thus producing the telegraph 
noise. Therefore, the noise spectral density S nn (ui) at 
frequencies lj <C <C Q, is 



S n n(<jj) 



T 



(A/„) 2 /2r s , (St n 



1 + (UJ/2TS) 2 



+ -jrSn, (3) 



where the term {5t n /T)S n is due to the shot noise. 
It is easy to see that at low frequency, ui <C Tg, 
the ratio of the telegraph and shot noise contributions 
(5t n /T)(AI n ) 2 /2T s S n ~ l/t col T s is always large in our 
case. 

For the phase shift p w ir the QND regime is still 
realized, and therefore Eq. (|3|) for each detector noise is 
still valid. However, since the detectors now measure the 
opposite qubit states, the cross-correlation noise changes 
sign, 



St a St b AI a AI b /2T s 

t 2 i + (w/2r s ) 2 ■ 



(4) 



where 



sign is for p fa and 



' is for ip sa 

7T. The noise correlation factor S a b{0)/\/ S aa (0)S bb (0) 
is obviously close to ±1, describing almost full correla- 
tion/anticorrelation, when the shot noise term in Eq. (|3]) 
is much smaller than the telegraph noise. 

To find the switching rate Ts, we calculate the "prop- 
agators" p^' to (t + t) in Eq. @, essentially rederiving 
Eqs. Q and (|4|) in the fully quantum way. Because of 
assumed weak coupling (7T -C M a , b -C 1) these den- 
sity matrices at time t = to + r can be represented 
in the Bloch coordinates as z — A(t)cos[ftt — 4>(t)], 
y = Tr(<jyp) = A(t) sin[fJt — <fi(t)] with slowly changing 
amplitude A(t) and phase 4>(t): 

A = -(A/T)[M a sm 2 (-cj)) + M b sin 2 (p - <j>)\ - 7 A/2,(5) 
4> = (l/2T)[M a sin(-20) + M b sm(2(p. - 2(f)], (6) 

where we assumed 5t ab <C T, so that periodic instan- 
taneous dephasings of magnitudes M a and M b happen 
at the phases —<$> and p — (j>. Assume now to = [so 
that f a (t ) = 1] and choose \i) = |1); then A(0) = 1 
and (f)(0) = 0. For \<p\ < 1 the solution of Eq. © is 
simple, and 4> saturates at 4> st = pM b /(M a + M b ) expo- 
nentially with time constant t co \. This value can be in- 
serted into Eq. (O because evolution of the amplitude A 
due to measurement is much slower than t~, , resulting in 
A(t) = exp[-(t/T)p 2 M a M b /(M a +M b )- 1 t/2}. It is easy 
to see that the evolution starting with the state |2) leads 
to 4> shifted by 7r, but the same A(t), and the same contri- 
bution into K zz (t, 0) in Eq. @. Calculating now S aa (oj) 
via K zz (t, 0), and using approximation | cos0 st | « 1 be- 
cause \4> s t\ < M -C 1, we obtain the formula, coinciding 
withEq. © withTs = {l/2T)p 2 M a M b /(M a +M b )+~//4. 
Calculation of S bb (uj) is fully similar, while to obtain 
S ab (oj) in the form ^ with the same switching rate Ts 
we also need approximation | cos(0 — 4> s t) \ ~ 1- 
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FIG. 2: (Color online) Numerical (solid lines) and analyti- 
cal (dashed lines) dependence of zero- frequency detector noise 
S aa (0) and cross-noise S a b(0) on the phase shift ip between the 
bias voltage combs for several values of the pulse width 5t a< h, 
and also for the harmonic biasing. Almost complete noise an- 
ticorrelation at p — ±7r indicate persistent Rabi oscillations. 



To account for small non-zero pulse widths St a ^ b , 
we can still use Eq. ©, but averaging over (j) within 
the pulse widths in Eq. ([5]) leads to the extra factor 
cxp[-(27r) 2 (M a <5i 2 + M b St 2 )t/12T 3 } in A(t) and corre- 
sponding increase of T$- A small frequency mismatch 
A£l = O — FIr would lead to the extra term Afl in Eq. 
©, so that (j) st becomes shifted by An T/(M a + M b ). In 
the case when the shift between the measurement combs 
is close to the half-period, p> should be obviously replaced 
by p ± 7T. Taking into account these changes, we repro- 
duce Eqs. ([3]) and (|4|) with the switching rate 



T s = 



<p 2 M a M b + (Af2T) 2 M a St 2 a 



2T(M a + M b ) 



6T 3 /tt 



M b 6t 2 _ + 1 
4 



for (p -C 1, where (p — min(|(p|, \p> ± 7r|). Assuming com- 
parable measurement parameters for both detectors, we 
see that the telegraph noise at zero frequency greatly ex- 
ceeds the shot noise contribution in Eq. ([3| if (p <C 1, 
5t a: b < T, |AO|T < M a>b , and 7 T < M a>b . This is the 
condition for the validity of our analytical results. 

Figure 2 shows zero- frequency spectral densities S aa (0) 
and S a b{0) as functions of the phase shift p for several 
values of the pulse width 5t a .b, assuming negligible Af2 
and 7. Solid lines are the numerical results calculated 
via Eq. ([2|), while the dashed lines show analytical re- 
sults using Eqs. ©, an d Q). Overall the analytics 
is very close to the numerical results (almost coinciding), 
except for S^O) near p — ±7r/2, where the analytics is 
discontinuous because of the sign change in Qj. We nor- 
malize the noise S' aa (0) by the shot noise S a of constantly 
biased detector, so that the shot noise contribution in 
this normalization is 5t a /T. Similarly, the cross-noise 
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Sab(O) is normalized by \/ S a Sb- The numerical results 
in Fig. 2 confirm almost full correlation of the detector 
noises at ip ~ and almost full anticorrelation at ip ~ ±7r 
[lfij . The peaks become higher and narrower for shorter 
pulse durations St a ^- The results in the used normaliza- 
tion are practically insensitive to the qubit-detector cou- 
pling 7™f as (assuming M a ^ <C 1). Non-zero detuning Aft 
and/or extra dephasing 7 make the peaks in Fig. 2 lower, 
while not affecting their width; this lowering is less sig- 
nificant for stronger coupling 7™. We have also checked 
numerically that the frequency dependence of the noises 
at ui <C ft is close to the analytical results ([3]) and ((4]); 
extra peaks as well as significant imaginary component 
of S a b(uj) appear at uj « ft and overtones of ft. 

Obviously, the analysis and results change only triv- 
ially if flu/ ft is an integer or close to an integer. In a 
real experiment with QPC detectors the best measure- 
ment mode is to apply two bias voltage pulses with op- 
posite polarity per Rabi period for each detector. In this 
case the average bias voltage is zero that helps to keep 
zero bias between the pulses. The average current in each 
detector is then also zero, simplifying the noise measure- 
ment. For such mode St a ,b in Eqs. ((3]) and (|4]) should 
be replaced by 25t a: b, while in Eq. ([7]) the measurement 
strengths M a ,b should be doubled (no change for 5t a ,b)- 

Now let us discuss why experimental observation of the 
noise dependence of Fig. 2 would indicate persistent Rabi 
oscillations. Correlation of the noises for <p could be 
alternatively explained by the qubit localization in either 
state |1) or |2). However, the anticorrelation for ip w ±ir 
is possible only if the qubit oscillates persistently. More- 
over, these oscillations should be synchronized with the 
measurement combs, because for persistent Rabi oscilla- 
tions with a random phase one would expect dependence 
S a b(0) oc cost/?, which is very different from the peaked 
dependence in Fig. 2. One may also worry that the noise 
dependence of Fig. 2 could be alternatively explained by 
the driven Rabi oscillations (between the energy eigen- 
states) caused by presence of a voltage with resonant fre- 
quency. However, both energy eigenstates produce no 
signal in the detectors; therefore the driven Rabi oscil- 
lations could only reduce the discussed noise correlation 
and cannot be used for an alternative explanation (notice 
also that both stable phases of the persistent oscillations 
are insensitive to the microwave drive oc cos fljit) . Un- 
fortunately, measurement of only zero-frequency noise is 
insufficient to demonstrate ^100% amplitude of the per- 
sistent Rabi oscillations (observed in [3j). However, if the 
switching rate T$ can be measured either in frequency or 
time domain, then by comparing S a b(0) with Eq. Q one 
can check the amplitude of synchronized oscillations. 

If the stroboscopic biasing is replaced by harmonic bi- 
asing [18]: f a = cos(Sli), f b = cos(flt - ip), then S ab (0) 
still depends on the phase shift ip (see Fig. 2); however, 
there are no more peaks and the noise magnitude is rela- 
tively small. The numerical results at weak coupling can 



be fitted as S ab {0) = 1. 18 AI a AI b cos ip/{Ya Cas + 7 £ ieas ) 
(actually, the (^-dependence is slightly more peaked than 

COSly9). 

In our analysis we have neglected the noise from am- 
plifiers, which can be simply added and is not expected 
to depend on ip. The main effect of the neglected thermal 
noise in the QPCs is a small contribution to the dephas- 
ing 7 between the pulses. A weak energy relaxation in 
the qubit can also be easily taken into account. 

For numerical estimates let us assume QPCs with 
I a b — 100 nA, AI a ^b/I a ,b — 0.1, symmetric biasing 
with 6t at b/T ~ 0.1, and Rabi frequency flu/2ir ~ 2 
GHz. Then the collapse ("attraction") time t co \ ~ 2 
ns is few Rabi periods, while the switching rate is T$ — 
(^ 2 /15ns + 1/120 ns+(Aft/ft) 2 /6ps + 7/4. Therefore we 
need the dephasing time Ti = 1/7 to be longer than only 
few ns to have significant correlated telegraph noise, and 
its ratio to the shot noise contribution for (p = Aft = 
is crudely min(60, T2/0.5ns) (5 times smaller for the nor- 
malization of Fig. 2). These figures show that the ex- 
periment is doable using the present-day semiconductor 
technology 0, [l9| . The experiment can also be realized 
with the superconducting qubit setup of Ref . Q . In com- 
parison with the experiment [3J it would also demonstrate 
partial synchronization of the persistent Rabi oscillations 
without use of a much more complicated quantum feed- 
back. 
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